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SPECTRAL DYNAMICS OF THE VELOCITY GRADIENT FIELD IN 

RESTRICTED FLOWS 

HAILIANG LIU AND EITAN TADMOR 


Abstract. We study the velocity gradients of the fundamental Eulerian equation, dtu+ 
u -Vu = F, which shows up in different contexts dictated by the different modeling of 
F’s. To this end we utilize a basic description for the spectral dynamics of Vu, expressed 
in terms of the (possibly complex) eigenvalues, A = A(Vm), which are governed by the 
Ricatti-like equation A* + u • VA + A^ = {l,VFr). 

We focus our investigation on four prototype models associated with different forcing 
F, ranging from simple linear damping and a viscous dusty medium models to the main 
thrust of the paper — the restricted models of Euler/Navier-Stokes equations and Euler- 
Poisson equations. 

In particular, we address the question of the time regularity for these models, that is, 
whether they admit a finite time breakdown, a global smooth solution, or an intermedi¬ 
ate scenario of critical threshold phenomena where global regularity depends on initial 
configurations. 

Using the spectral dynamics as our essential tool in these investigations, we obtain 
a simple form of a critical threshold for the linear damping model and we identify the 
2D vanishing viscosity limit for the viscous irrotational dusty medium model. Moreover, 
for the n-dimensional restricted Euler equations we obtain [n/2] -|- 1 global invariants, 
interesting for their own sake, which enable us to precisely characterize the local topology 
at breakdown time, extending previous studies in the n = 3-dimensional case. Einally, 
as a forth model we introduce the n-dimensional restricted Euler-Poisson (REP)system, 
identifying a set of [n/2] global invariants, which in turn yield (i) sufficient conditions for 
finite time breakdown, and (ii) characterization of a large class of 2-dimensional initial 
configurations leading to global smooth solutions. Consequently, the 2D restricted Euler- 
Poisson equations are shown to admit a critical threshold. 
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1. Introduction 

It is well known that the velocity gradients in a turbulent flow are larger than their mean 
gradients by at least a factor of order with Rs being the Reynolds number based 

on internal length and velocity scales. Fluctuation gradients are limited by the mean 
flow and contribute a dominant portion of the kinetic energy dissipation, but otherwise 
they contribute nothing to the mean transport of momentum because of the linearity of 
the viscous stress term in the Navier-Stokes equations. Consequently, much research has 
been directed at gaining a better understanding of the velocity gradient held, Vm, which 
is completely dictated by the vorticity in incompressible flows. 

Motivated by such questions, the goal of this work is to present new observations on 
the velocity gradients for a general class of so called restricted flows, where the velocity 
held, u, is governed by the Newtonian law, 

(1.1) dtU + u-Vu = F, a; e IR”", 

with F being a general forcing acting on the flow. Different regimes of the flow are 
modeled by different F's. A key issue in this line of research is the control of the velocity 
gradient Vm, and a classical approach in this context, is to consider linear combinations 
of the entries of Vm, controlling physically relevant quantities like vorticity, divergence, 
etc, see i, 0, 0. 

The novelty of the analysis taken in the present article is the use of the eigenvalues of the 
velocity gradient held. The eigenvalues, A = A(Vn), exhibit of course a strong nonlinear 
dependence on the entries of Vu, and are shown to play a crucial role in governing the 
behavior of the flow. Indeed, the dynamics of these eigenvalues A(Vm), is shown, in §3, 
to be governed by Ricatti-like equation 

dtX + u-VX + A^ =< I, VFr >, 

with l{r) being the left (right) eigenvectors of Vm. Equipped with this description for the 
spectral dynamics of Vm, we turn to study several physical models with different forcing, 
outlined in §2 and analyzed in §4-§7. 

We focus on four prototype models in this paper. The first two are a simple linear 
damping model studied in §4, and a viscous dusty medium model in §5. Next, the main 
thrust of the paper is devoted to the restricted models of Euler equations in §6, and 
in §7 we introduce the so called restricted Euler-Poisson equations as our fourth model 
problem. We focus our attention on the question of time regularity for these models, 
that is, whether they admit a finite time breakdown, a global smooth solution or an 
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intermediate scenario of critical threshold phenomena where global regularity depends on 
the initial conhgurations as in e.g. 0. 

The question of time regularity is of fundamental importance from both mathematical 
and physical points of view, and a considerable effort is still being devoted to this issue for 
both compressible and incompressible Euler equations. Consult |^, 10, n, i I, n, n 
for a partial list of recent references. In particular, the possible phenomena of hnite 
time breakdown for 3D incompressible flows signihes the onset of turbulence in higher 
Reynolds number flows. Several simplihed models for 3-D Euler equations were proposed 
to understand this phenomena, see for a restricted dynamics model, for a vorticity 
dynamics model, for a so called distorted model as well as a stochastic model in 0. 


The paper is organized as follows. 

After introducing the basic spectral dynamics Lemma 3.1 in §3, we begin our discussion 
of the time regularity with the simple linear damping forcing model in §4. 

In §5 we deal with viscosity forces, where we study the irrotational viscous dusty 
medium model, and identify its 2D vanishing viscosity limit. Here, the spectral dynamics 
offers us a novel approach at the level of the velocity held, m = V0 — an alternative to 
the classical notion of viscosity solutions for Hamilton-Jacobi equations at the level of the 
potential 0. Spectral dynamics serves as an essential tool in our approach, most notably 
the use of a key a priori L^-contraction estimate expressed in terms of the unintuitive 
nonlinear quantity A 2 — Ai, with Aj = Xi{'Vu),i = 1,2 being the two real eigenvalues of 
Vm. 

In §6 we use the spectral dynamics to revisit the restricted Euler models introduced 
by Vieillefosse in [0. The so called restricted Euler equations (RE for short) refer to 
a localized model of the Euler/Navier-Stokes equations, where the usual global pressure 
forces are replaced by their local, isotropic trace. We study the time regularity of the 
general n-dimensional RE equations, extending the previous studies in [04|, for the 
special n = 3 case. Here, we enjoy the advantage of using the spectral dynamics of 
such n-dimensional hows, which enables us to identify a large set of at least [n/2] -|- 1 
independent integrals of the motion. Using these [?7,/2]-|-l global invariants, interesting for 
their own sake, we precisely characterize the hnite time breakdown for the n-dimensional 
RE equations. 

We note in passing that the RE model has been an appealing candidate for describing 
the dynamics of the local velocity gradient, BB- Despite this restricted approximation 
to the pressure, the RE equations can still describe the local topology of Euler equations 
and capture certain statistical features of the physical how. In this spirit we introduce, in 
§7, a restricted model for the Euler-Poisson system, so called REP equations. For general 
n-dimensional REP equations we obtain a set of at least [n/2] global invariants, which in 
turn yields 

(1) Sufficient conditions for the hnite time breakdown in n-dimensional REP equations. 
Moreover, we characterize the precise local topology of the how at breakdown time; 

(2) Sufficient conditions for a large class of 2D initial conhgurations leading to the 
existence of global smooth solution for 2D REP equations. 

We point out that though RE model was sometimes argued for its unphysical hnite 
time singularity, our REP model does support the global smooth solutions. In particular, 
it follows that the 2D REP equations admit a critical threshold, distinguished between 
initial conhgurations leading to either the hnite time breakdown or the global smooth 
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solutions. We refer to [Q for a detailed study of such phenomena for this 2D REP 
model. 

In §8 we discuss possible extensions of the results obtained in this work and we comment 
on some remaining open issues. Finally, in the Appendix we revisit the spectral dynamics 
of the general nxn RE models from yet another perspective of a trace dynamics, extending 
the study of traces, tr(Vu)'^, k = 1, ■ ■ ■ ,n, initiated in for the special n = 3 case. 


2. Basic equations — four prototype models 


In what follows we shall require the equations governing the dynamics of a fluid in both 
the Eulerian and Lagrangian forms. We shall study the flow of a fluid which initially at 
t = 0 occupies the whole space IR"^ for arbitrary dimension n, although only the cases 
n = 2 and n = 3 have a clear physical meaning. 

Let a Cartesian coordinate system be fixed in IR"'. We denote by a the initial position 
of a fluid particle. The motion of the fluid is assumed to be given, if for any a G IR" the 
position x{a, t) G IR” of the fluid particle is known for all a G IR” and for all t G IR^, with 
x{a,0) = a. Further, a —> -^x{a,t) = u{x,t) is a Lagrangian velocity field at the time t. 
The Lagrangian equations of the dynamics of a fluid amount to 


(f 


F, 


where F denotes the forcing acting on the fluid. The corresponding Eulerian equations 
in the standard form read 


( 2 . 1 ) 


dtu + u ■ Vm = F, a; G IR", 


where ^ ^ + u ■ V is the Lagrangian derivative. Equation (p2.1|) shows up in a variety 

of contexts dictated by the different modeling of F’s. 

Differentiation of the above equation with respect to x gives the relation for the local 
velocity gradient tensor, M := Vu, 

(2.2) dtM + {u ■ V)M + M^ = VF. 


The central issue of interest here is to control the local velocity gradient tensor in (|2.2|) 
and to clarify whether the associated distortion matrix, T := dx/da, remains nonsingular 
as time evolves. In particular, one is interested to know whether there is a finite time 
breakdown, a global smooth solution or an intermediate scenario of critical threshold 
phenomena with a conditional breakdown, consult 


e.g. 


15 


In the remainder of this 
section we discuss four prototype examples associated with different forcing F. In the 
following sections, we will follow the spectral dynamics of the velocity gradient tensor 
associated with these four examples to demonstrate the different phenomena of global 
regularity, finite time breakdown as well as the existence of critical threshold. 


2.1. Linear damping. Consider a model of the form 

(2.3) dtU + u ■ Vu = Cu, X G IR”. 

Here we deal with the simple forcing, F = Cu, where C is a constant matrix. The 
corresponding local velocity tensor M = Vu solves 

(2.4) dtM + {u ■ V)M + = CM. 
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In we use the spectral dynamics of M to show that there exists a critical threshold 
depending on the choice of the matrix C. 


2.2. Irrotational viscous flow. Next we consider the viscous forces, F := i/Am, which 
leads to the so called viscous dusty medium model 


(2.5) 


dtu + u ■ Vm = i/An, X G IR”, 


where z/ > 0 is a viscosity amplitude. Other suggested names are Burgers system 0, 
Hopf system, Riemann equation (for n = 1). Zeldovitch proposed to consider this 
system as a model describing the evolution of the rarehed gas of non-interacting particles. 

The Hopf-Cole transformation, u = —2v'^\Log{'ipy\, links the Burgers system to the heat 
equation 

dt'ip = 


provided the initial data, mq = m(x, 0), admits the form uq = —2z/V[/o5f(V’o)] for some 
positive '^0 = 0) (this is available for n = 1). The corresponding local velocity 

gradient held satishes 


( 2 . 6 ) 


8,M + (u ■ V)M + Nf = i/AM, x £ IR”. 


We focus our attention on solutions to the 2D irrotational viscous flows, u = and we 
use the spectral dynamics of M in order to study the inviscid limit, u = In 

particular, the limiting 2D irrotational flow is shown to be a weak solution of 


dtU + u ■ Vu = 0, u =-. V0, 

which is interpreted through the Eikonal equation dt(j) + | V0P/2 = 0. 


2.3. Restricted Euler/Navier Stokes equations. For the forcing involving viscosity 
and pressure, we consider the Navier-Stokes equations of incompressible fluid flow in n 
space dimensions, which can be expressed as the system of n -|- 1 equations, 

(2.7) dtu -|- u ■ Vu = uAu — Vp, x G IR"", t > 0, 

(2.8) V-M = 0, 

(2.9) u{x,0) = uo{x). 

Here u is the fluid velocity, p is the scalar pressure, and z/ > 0 is the reciprocal of the 
Reynolds number. When the coefficient z/ vanishes in (|2.7|) , we have the incompressible 
Euler equations. Here we only discuss fluid flows occupying the whole space so that the 
important effects of boundary layers are ignored. In most applications, u is an extremely 
small quantity, typically ranging from 10“^ to 10“® in turbulent flows. Thus one can 
anticipate that the behavior of inviscid solutions of the Euler equations with z/ = 0 is 
rather important in describing solutions of the Navier-Stokes equations when z/ is small. 
The local velocity gradient tensor solves 

(2.10) dtM + (u ■ V)M + = z/AM - (V 0 V)p. 

Taking the trace of M and noting trM = V ■ m = 0 one has 

(2.11) trM'^ = —Ap. 
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This gives p = —A ^{trM'^). The second term in ( |2.1CI| ) therefore amonnts to the n x n 
time-dependent matrix 

(V ® V)A-\trM^) = R[trM% 

Here R[w] denotes the nxn matrix whose entries are given by {R[w])ij 
Rj denote the Risez transforms, Rj = i.e., 


:= RiRj{w) where 




for 1 < j < n. 


This yields the eqnivalent formnlation of NS eqnations which reads 
(2.12) dtM+{u-V)M+ M‘^ = vAM+ R[trM‘^] 

snbject to the trace-free initial data 

M(-, 0) = Mo, trMo = 0. 

Note that the invariance of incompressibility is already taken into account in ( p.l2|) since 
dttrM = 0 and hence trM = trMg = 0. The inviscid case z/ = 0 in (|2.12|) gives the 
corresponding Euler equation. It is the global term in the above equations, i?[frM^], 
which makes the problem rather intricate to solve, both analytically and numerically. 
Various simplifications to this pressure Hessian were sought, see, e.g. m 0 ini- 
Here we focus our attention on the so called restricted Euler equations proposed in 
as a localized alternative to (p.l2|) . 

Specifically, we consider a gradient flow, M, governed by 




(2.13) 


dtM -|- {u ■ V)M -|- M^ = —trM'^In 


n 


We observe that as in the global model, the incompressibility is still maintained in this 
localized model, since trM‘^ = tr[trM‘^Inxn/'n] implies dttrM = 0. 

For arbitrary n > 3, we use the spectral dynamics of M in order to show a finite time 
breakdown of ( p.l3D , generalizing the previous result of 13^. The finite time breakdown 


follows in once we identify a set of [n/2] -|-1 global invariants in terms of the eigenvalues 
of M. Moreover, the precise topology of the flow at the breakdown time is studied in §|^. 
Finally, in the Appendix we study the spectral dynamics of the general nxn problem from 
yet another perspective, extending the study of traces, trM^, k = 1, - ■ ■ ,n initiated in 
H for the special n = 3 case. 


2.4. Restricted Euler-Poisson equations. The Euler-Poisson equations 


(2.14) 

Pt + V ■ {pu) = 0, X E IR"", t E IR^ 

(2.15) 

{pu)t -h V ■ {pu ® u) = kpVcj), 

(2.16) 

A0 = p, X E IR"', 


are the usual statements of the conservation of mass, Newton’s second law, and the Poisson 
equation defining, say, the electric field in terms of the charge. Here fc is a scaled physical 
constant, which signifies the property of the underlying forcing, repulsive if fc > 0 and 
attractive if fc < 0. The unknowns are the local density p = p{x,t), the velocity field 
u = u{x, t), and the potential 0 = (j){x, t). 

If follows that 


dtu + u ■ Vu = 
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where the forcing F = kVcf) is the gradient of potential governed by the Poisson equation 
( p.l6| ). Differentiation yields a local velocity gradient tensor which solves 

dtM + u ■ VM + = k{y ® V)0 = kR[pl 

where the coupling enters through the global term fc-R[p], with density p governed by 

dtp + M ■ Vp + ptrM = 0. 

Passing to Lagrangian coordinates, that is, using the change of variables a i—> x{a, t) with 
x{a,t) solving 


dx 

dt 


= u{x, t), x{a, 0) = a, 


then Euler-Poisson equations become the coupled system 


(2.17) 


M + M^ = kR[p], 


d 

dt 


:= dt + u-V, 


—p + ptrM = 0, 
dt 


(2.18) 

subject to initial condition 

(M, p)(-, 0) = (Mo, po). 

Again, it is the nonlocal term, R[p], which is the main obstacle, in the multi-dimensional 


setting n > 1, in the investigation of the Euler-Poisson system, see e.g. [21 


In this paper we introduce the corresponding restricted Euler-Poisson system 


dtM + u ■ VM + M^ = -pk 

n 

dtpu -Vp + ptrM = 0, 


nXTL’) 


(2.19) 

( 2 . 20 ) 

subject to initial data 

(M, p)(-, 0) = (Mo, Po)- 

In we use the spectral dynamics of M in order to study the time regularity for this 
restricted Euler-Poisson model. Here we give a sufficient condition for the global existence 
of the 2D solutions which applies, for example, for a class of initial conhgurations with 
sufficiently large vorticity |a;o| >> 1. With other initial conhgurations, however, the hnite 
time breakdown of solutions may - and actually does occur. Hence global regularity 
depends on whether the initial conhguration crosses an intrinsic, 0(1) critical threshold, 
and we refer to for a detailed study of the 2D critical threshold phenomena in this case. 
Moreover, for arbitrary n > 3 we obtain a family of [n/2] global invariants, interesting for 
their own sake, with which the local topology of hnite time breakdown is also characterized 

in §0. 


3. Spectral dynamics of the Velocity Gradient Field 

Let us rewrite the basic equation of velocity gradient held M as 

(3.1) a^M + (m • V)M + M^ = VF, 

where VF is a matrix involving spatial derivatives of the forcing. 

It is usually difficult to quantify directly all entries in the velocity gradient tensor, 
M, and instead, suitable linear combinations like divergence, vorticity play a distinctive 
role in analysis. Here we show the special role played by the eigenvalues of the velocity 







H. Liu and E. Tadmor 


gradient tensor, A(M), in governing the entries of M, and we note in passing, the strong 
nonlinear dependence of A(M) on the entries of M. Consnlt, for example, the nonintnitive 
L^-contraction for 2D dnsty medinm model derived in (|5.9|) below. 

The following lemma is in the heart of matter. 


Lemma 3.1 (Spectral dynamics). Consider the general dynamical system B) associ¬ 
ated with arbitrary velocity field u and forcing F. Let A(M) he a (possibly complex) eigen¬ 
value of M associated with corresponding left(right) eigenvector I (r). Then the dynamics 
of X{M) is governed by the corresponding Ricatti-like eguation 


dt\ + u-V\ + A^ =< /, VFr > . 


Remark 3.2. If F = 0 one has the same eqnation for A as for M with time-independent 
eigenvectors, thus M{t) are isospectral. The difficulty lies in the eigenstructure induced 
by the forcing < /, VFr >. 

Proof. Let the left (right) eigenvectors of M associated with A be l{r), normalized so that 
Ir = 1. Then one has 


Mr = Ar, IM = XL 

Differentiation of the hrst relation with respect to t gives 

dtMr Mdtr = dtXr + Xdtr. 
Multiply I on the left of the above equation to obtain 

IdfMr Xldtr = dtX -\- Xldtr, 


hence 

IdtMr = dtX. 

Similarly differentiation of the relation Mr = Ar with respect to Xj leads to 

djMr -|- Mdjr = djXr -|- Xdjr. 

Multiply on the left by luj with /r = 1 to get 

lujdjMr = UjldjXr = UjdjX. 

Therefore 

lu ■ S/Mr = u ■ VA. 

A combination of the above facts together with lM‘^r = A^ gives 

dtX + u-VX + A^ =< /,VFr > . 


This completes the proof. 


□ 
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4. Critical Thresholds for Linear Damping 
Consider the convective eqnation 

dtu + u ■ Vm = Cm, u{x, 0) = Uo(x), 

with a simple linear damping C being a constant matrix. The gradient tensor satishes 

dtM + u ■ VM + = CM, 

which snggests that the eigenvalnes solve 
(4.1) ^tA + M-VA + A^ = cA, 

where c(t) = CM{t) '■= ICr. Here l{r) are the left (right) eigenvectors of M associated with 
the eigenvalne A. Along the particle path x = x{a,t), dehned by 

■^x(a,t) = u(t, x(a,t)), x(a,0) = a, a G IR"', 

at 


the Ricatti-type A-eqnation amonnts to 

d 


dt 


A + A^ = c(f)A. 


The solntion can be expressed in terms of c(t) as 

mm 


b{t) := exp ( / c{T)dT . 


X(t) = - . , 

1 + A(0) /o 6(r)dr 

^From the above formnla it follows that 

Lemma 4.1. Consider the eigenvalue equation (M) with initial data A(0). 

(1) If Im{X{0)b{t)) 7 ^ 0, then its solution remains regular for all time; 

(2) If Im{X{0)b(t)) = 0, then its solution remains bounded as long as 


(4.2) 


Re ( A(0) / b{T)dT ) > —1. 


For the simple example of a scalar damping, C = —fllnxn, /9 > 0, one has 


dX 

dt 


+ X^ = -/3A, 


with a solntion (corresponding to b{t) = e ^*) given by 

A(a,0)e-^* 


X{t) = 


l + A(a,0)/9 ^(1 — e ht'j 


This solntion is bonnded from below for all time if and only if A(q!, 0) is either complex 
or 

inf A(a, 0) > —ft, 

which is a very simple form of a critical threshold phenomena. For more general examples 
of critical threshold phenomena, consnlt [l^, |l^, ^ and the stndy in §0 below. 
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5. Irrotational viscous flow 

In this section we deal with viscous dusty medium flow, u := u'', governed by 

(5.1) dtu + u-Vu = u{x,Q) = uq{x). 

The velocity gradient tensor M :='Vu satisfles 

(5.2) dtM + u ■ VM + = uAM. 

It follows that if the initial velocity is irrotational, V x rq = 0, then the flow remains 
irrotational, V x m = 0. 

Lemma 5.1 (Viscous Spectral Dynamics). Assume that the flow is irrotational V xuq = 
0. Then the real eigenvalues A = A(Vr) satisfy 

dtX + u ■ VA + A^ = z/AA + Q. 

Here Q satisfies the constraint 

a{Xm,in -X) <Q < a{Xmax “ A), A := I I A(Vm), 

^ min ^ 

where a is given by 

a := 2u ^ dkV^dkV > 0 

k 

and r being the right eigenvector ofVu associated with X. 

Proof. Let l{r) be the normalized left (right) eigenvectors of M associated with the eigen¬ 
value A, then one has 

dtX + u-VX + X^ = ulAMr. 

Observe that M is symmetric due to the fact that V x w = 0, and consequently A are all 
real quantities. Differentiation of IM = XI with respect to x twice gives 

AIM + 2V/ ■ VM + I AM = AXl + 2VA ■ V/ + AA/, 

which upon multiplication against r on the right leads to 

lAMr = AA + 2 [(VA ■ V/)r - (V/ ■ VM)r]. 

Here the differentiation operators apply component wise, e.g., V/ ■ VM = ^flOkM. 
On the other hand it follows from Mr = Ar that 


VMr = VAr -|- AVr — MVr. 

This gives 

(V/ ■ VM)r = V/ ■ VAr + XVI -Vr-Vl- MVr. 
A combination of the above facts yields 


Q = 2u -A dkldkr + dfiMOkr . 

^=1 k=\ 

Since the flow is irrotational we have M"*" = M and I = r~^, with upper-index denoting 
the transpose. The second term in Q is then bounded by 

n n n 

^min E dkr^dkr < ^ dflMdkr < X 

max E dkr^dkr, 

k=l k=l k=l 
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which completes the proof. □ 

Here the question of interest for us is the convergence of u = as z/ —>■ 0. To answer 
this question, it suffices to show the precompactness of the family {u^}u>o- It is here 
that we take advantage of the spectral dynamics of the velocity gradient tensor, Vu. For 
the 2D case we shall show the precompactness via several lemmata. We start with the 
essential 


Lemma 5.2 (L^-Contraction). Let A,, i = 1,2, be two (real) eigenvalues of velocity gra¬ 
dient field in If (A 2 — Ai)(0) G L^(IR^), then 

11(^2 - Al)(t)||il(|R2) < ||(A2 - Ai)(0)||^1(|r2). 

Proof. In the 2D case we have A^m = Ai < A 2 = \max- Setting rj = X 2 — Xi one has from 
Lemma 4.1 

dtT] + u ■ Wr] + r]{Xi + A2) < vlX-q. 

Observe that V ■ n = Ai + A 2 which yields 

dth + V ■ (qu) < uArj. 

Spatial integration gives the estimate for 77 = A 2 — Ai > 0 as asserted. □ 


Next, from Lemma 5.1 we see that the largest eigenvalue Ai, 
inequality 

^zAmax T If ■ VAjiiax T -^max — f^^'^max 
and by a comparison principle we obtain 


satisfies the differential 


(6.3) 


Amax(t) ^ 


1 ^1 
A(0)-1 + f “ T 


We note in passing that this, combined with the symmetry of gradient field Vn, is equiv¬ 
alent to the one-sided entropy-type estimate supn^n^]^ < 1/t, which coincides with 

the well known semi concavity property in the context of convex Hamilton-Jacobi equa¬ 


tions, see, e.g., p5 


The above one-sided bounds enable us to establish the following. 


Lemma 5.3 {BV Bound). Consider the dusty medium equation (EIP with compactly 
supported initial data Uq = 0 ) such that i|no||By(iR 2 ) is bounded uniformly in v. 

Then the corresponding velocity, u’^, satisfies 

\\u’'{-C)\\BV(\Rfi < Const. 

Moreover, for ti,t 2 >ll we also have 

( 5 . 4 ) ||M'^(a:, ^2) - ti)||Li(iR2) < Const.|t2 - 

Proof. The one-sided upper bound for Xmax, implies that the positive part of the 

divergence, + ny)+ = (Ai -|- A 2 )+ is bounded. We observe that Ai, A 2 are essentially 
supported on a hnite domain in the sense of their exponential decay outside a finite region 
of propagation, and hence /|R 2 ('n 3 ; + Vy)+ < Const. This, combined with J|p 2 + Vy = 0, 
yields that Ai -|- A 2 = -1- G L^(IR^). Augmented with the fact that A 2 — Ai G L^(IR^) 

we conclude 


(5.5) 


A.gL1(IR2), 7 = 1,2. 
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This gives 


\^Vu^\^dxdy = / \^diag{\i,\2)\\dxdy < oo, 


'ir2 




with the usual matrix norm, || ■ ||, defined as ||M|| =: sup||^||=]^ \^i\- since M = Vu 

is symmetric, there exists a unitary matrix U such that U~^MU = diag{Xi, A 2 ), and hence 

\\\/u\\ = \\U^MU\\ = \\dtag{X,,X 2 )\\. 


Thus, the BV bound of follows from ( b-5|) . To estimate the modulus of continuity in 
time, we multiply equation ( b-l|) by a smooth test function ^lJ G and use the spatial 
BV bound to obtain 


'iIj{x){u{x, ^ 2 ) — u{x, ti))dx 


'IR^ 


< Const.(f 2 - ^i)(|'0|oo + |A'0|). 


This inequality and the BV estimate combined with Kruzkov’s interpolation theorem [^, 
page 233] yield □ 


In order to identify the vanishing viscosity limit, limj^^ow^', we introduce the notion of a 
weak solution for corresponding inviscid equation 

(5.6) dtu + u ■ Vu = 0. 

For irrotational flow, V x u = 0, one has u ■ Vu = V (|np/2), and the reduced inviscid 
equation (|5.6|) can be recast into the conservative form 

dtU + 'V [1 = 0 . 


The irrotational property of both viscous and inviscid flows suggests that there exists a 
potential 0 such that u = V0, where 0 solves the Hamilton-Jacobi equation 


(5.7) 


4>t + = 0, 0(x,O) = 0o- 


According to the classical theory of the Hamilton-Jacobi equation 0,i, there exists a 
unique continuous solution (j){x, t) to the above problem, expressed via the Hopf-Lax’s 
formula, [^, page 560], (l>{x,t) = miHygiR" {t\x — y\^/2 + 0o(l/)} • We make 

Definition 5.1. A measurable function u is called a weak solution of the inviscid equation 
( b.6|) if M = V0 with the potential 0 being the unique weak solution of the Eikonal equation 

(El- 

Equipped with this definition of a weak solution, we now turn to summarize our con¬ 
vergence results by stating 

Theorem 5.4 (Vanishing viscosity limit). Consider the dusty medium equation l \5.1\ ) with 
irrotational initial data m‘^(-, 0) G fl L°°(IR^) such that 

u''{x,0) ^ uo{x) in L^(IR^). 

Then, the local velocity u'^ converges to the unique weak solution of 457^ ), i.e., we have 
(5.8) u''{x,t) ^u{x,t) m L“([0,T];L^(IR2)), 

where u = Vcf) is the viscosity solution of the Eikonal equation (M)- 
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Proof. We begin by first assuming that u’^{x,0) is compactly supported in i?l/(IR^), uni¬ 
formly with respect to u. By Lemma 4.2, u'^ have uniformly bounded spatial variation, 

i.e., 

ll«''(•,^)llw(lR2) < Const. 

Hence t < T} is a bounded set in fl BV{\R^) and by Helly’s theorem it 

is therefore precompact in A'oc(IR^)- Note that is Holder continuous in 

time, and by Cantor diagonalization process of passing to further subsequence if necessary, 
(5.8) follows. This completes the corresponding proof for compactly supported BV initial 
data. The general case is justihed by standard cutoff and BV-regularization of arbitrary 
n L°°(IR^) initial data. 

It remains to show that the limit function u{x, t) satisfy the weak formulation. It follows 
from the equation for that 


dtu^ + V 


\u 


ul2 


= uAu’' 


We multiply this identity by fj{x,t) E C“(IR^) and integrate by parts to get 

,I2\2- 


' R? 


-'ipt'iA -Vip- 


\u 


dxdy = jy u'' A'tpdxdy. 


'IR^ 


Note that is embedded into L^(IR^) for the two-dimensional case. Thus passing to 


the limit i/ —0 one obtains the desired weak formulation. 


□ 


Remark 5.5. We would like to point out that the above convergence result can be obtained 
at the level of Hamilton-Jacobi equations. The equivalence between the weak entropy 
solutions to conservation laws and the viscous solutions to the corresponding Hamilton- 
Jacobi equations has been known in the literature, see e.g. 1,0,0,0 0. The point 
made here is that we obtain the compactness at the level of u by using the spectral 
dynamics of its velocity gradient tensor, Vm, which is independent of a maximum principle 
at the level of HJ equations. In particular, the 2D L^-contraction stated in Lemma 5.2, 
recast at the level of HJ equation (|5.7|), amounts to the nonintuitive apriori estimate 


(5.9) 




< 


Li 


\/ (^0o)^ J (0o) ) '^(0) • 


Ll 



6. Restricted Euler Dynamics 


6.1. Spectral dynamics and global invariants. We now turn to discuss the restricted 
Euler dynamics, which is a localized version of the full Euler/Navier-Stokes equation 

dtM + (m ■ V)M + M^ = pAM + R[trM^]. 

By the dehnition of the operator R, one has 


R[trM^] = V ® VA-^[trM^] 
where the kernel K{-) is given by 


V ® V 



K{x) 


yn\x\ 


2n 


(2-n)uj„\x\ 




y){trM^){y)dy, 


n=2, 
n > 2, 
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with uJn denoting the snrface area of the nnit sphere in n-dimensions. A direct compntation 
yields 


did.K * trM^ = ——h 


n 




+ 


k - y?^i3 - - yj) 


MR" 


Cd, 


\x - y\ 


n+2 


trM‘^{y)dy. 


This shows that the local part of the global term R\trM‘^] is trM"^/ninxn- We thus use 
this local term, trM^/ninxn, to approximate the pressure Hessian. The corresponding 
local gradient held then evolves according to the following restricted Euler model 

( 6 . 1 ) dtM + u ■ WM + M'^ = trM'^/ninxn- 


This is a matrix Ricatti equation which, as we shall see below, is responsible for the 
formation of singularities at hnite time. We note that with this local model, all particles 
evolve independently of each other. The mixing due to the global forcing in the general 
Euler dynamics, however, could prevent this type of hnite time breakdown. 

Nevertheless, as a local approximation of the pressure Hessian, the above model, so 
called restricted Euler dynamics, has caught great attention since hrst introduced in , 


because it can be used to understand the local topology of the Euler dynamics and to 
capture certain statistical features of the physical how. 

Consider a bounded, divergence-free, smooth vector held u : IR” x [0,T] —> IR”. Let 
X = x{a,t) denote an orbit associated to the how by 

dx 

— = u{x, t), 0 < t <T, x{a, 0) = a G IR". 

(JjL 

Then along this orbit, the velocity gradient tensor of the restricted Euler equations (|6.1|) 
satishes 

M' + M^ = trMynInxn, ' ■-= 

By the spectral dynamics lemma 3.1, the corresponding eigenvalues satisfy 


( 6 . 2 ) 


K + 


= 1 , 


, n. 


k=l 


This is a closed system which serves as a simple approximation for the evolution of the 
velocity gradient held. 

Let us start by revisiting the case n = 3, consult 


34, O, for which we will present below 


a complete phase-plane analysis expressed in terms of A's. Subtraction of two consecutive 
equations in ( p. 2 |) gives the following equivalent system 


[/?7.(Ai — A2)]^ -|- Ai -|- A 2 — 0, 
[/n(A 2 — A 3 )]' -l- A 2 -l- A 3 = 0, 
[/?7,(A3 — Ai)]' -|- A 3 -|- Ai = 0. 


Summation of these three equations and taking into account the incompressibility condi¬ 
tion, = Oi yields the following global invariant 

(Ai — A 2 )(A 2 — A 3 ) (A 3 — Ai) = Const. 

This invariant projected onto the phase plane (Ai, A 2 ), recast into 

(A 2 — Ai)(A 2 -I- 2 Ai)( 2 A 2 -I- Ai) = Const., 
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which serves as a global invariant of the 2x2 system 

(6.3) A( = [—+ 2 A 2 + 2AiA2]/3, 

(6.4) A 2 = [ 2 A 2 — A 2 + 2AiA2]/3. 

We then have three separatrixes passing throngh the origin, which is the only rest point 
in this case, 

Ai = A2, Ai = —2A2 ^ (A2 = A3), Ai = ~ 2 ^'^ ^ ~ 

The vector held in the phase plane is drawn in Fignre |6.1|. 


3-D Spectral dynamics 



Figure 6.1. 3-D Spectral dynamics of the Restricted Enler Eqnations 


Three special solntions corresponding to the separatrixes can be obtained explicitly. 
Consider, for example, the separatrix Ai = A 2 , for which Ai is necessarily a real solntion 
of the Ricatti eqnation 

a; = A?. 

The solntion, given by 


(6.5) 




Ai(q;, 0) 

1 — Ai(a, 0)f ’ 


is bonnded if and only if the real Ai is nonpositive, Ai(a, 0) < 0. 

Next if Aj(0) is complex then Aj remains complex later on. Let (Ai,A 2 ) be a complex 
pair with X 2 = Xi = ■y + Pi, then (|6.3|) , (|6.4|) recast into 


P' = -2Py, y = f + 


(6.6) 
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Solving the above 2 x 2 system gives the following invariant, (/9^ + = Const; it 

follows that the general solution passing the rest point ( 0 , 0 ) must be real, 


^ = 0, 7(t) 


7(0) 

1 -7(0)t’ 


which is reduced back to the hrst case in (^). Note that if the eigenvalues are complex. 


the Lagrangian trajectories are rotating, and if the eigenvalues are all real, the Lagrangian 
trajectories are just strain dominated, see Figure 


3—D rotational flow 



Figure 6.2. 3-D Rotational flow in Restricted Euler Equations 


We now summarize by stating the following 

Lemma 6.1 (Topology of flow in 3-D restricted Euler). A global invariant of the 3D re¬ 
stricted Euler equations l\6. is given by 

(6.7) (Ai — A 2 )(A 2 — A 3 )(A 3 — Ai) = Const. 

The three explicit solutions passing through origin are 

(A,. A.,, A3)(() = {(1,1, -2), (1, -2,1), (-2,1, 

1 — a{a, (Jjr 

All other solutions will develop finite time singularity. If an eigenvalue is complex, then 
the Lagrangian trajectories are rotating. 

Lemma EH deals with the 3-dimensional restricted Euler equations which were studied 
earlier in using a different approach based on trace dynamics, consult the Appendix 
below. Here we enjoy the advantage of being able to generalize our spectral dynamics 
approach taken in Lemma |6-1| to the arbitrary n-dimensional case. The global invariants 
in such n-dimensional systems are tied to a particular set of sequences of indices X = 
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{Xi,X 2 ,... } with each X^ being a sequence of pairs of different indices, {i ^ j), such that 
there exists an integer N := N{n) for which 


( 6 . 8 ) 




(Aj + Xj) — N ^ ^ Afc, WX's. 


k=l 


There are several ways of forming these pairs, (f, j), so that (| 6 . 8 |) holds. Here is one: 

• For even n = 2m we let = {i,j) = {{a{2k — 1), 2 ranging overall 

permutations a(-) so that (|6.8|) holds with iV(n)|{„ even} = 1; 

• For odd n = 2m + 1 we let X^.^ = (i,j) = {(o'(k) 7 ^ /(f(fc))}fc=i 2 n> ranging overall 
permutations a(-),ju(-) so that (| 6 . 8 |) holds with N(n)i^n odd} = 2 . 

The following lemma reveals the role such X-pairs play in forming global invariants for 
the restricted Euler system (|6.1|). 


Lemma 6.2. [Global Invariants] Consider the n-dimensional restricted Euler system ( \6.t\) 
subject to incompressible initial data, Ai( 0 ) = 0 . Then it admits the following global 
invariants in time: 

n 

(6.9) ^A,(f) = 0, 

i=l 

(6.10) (Aj(t) - Xj{t)) = Const. 


Proof. Summation of the equations in (|6.2|) over index i gives [X]r=i ~ 0’ which is 

combined with the incompressibility assumption Aj( 0 ) = 0 to yield X]r=i ~ 0 - 
For (16.101) we follow our previous argument in the 3D case. Subtracting the j-th equa¬ 
tion from i-th equation in ( b. 2 |) yields 

(6.11) T (A 2 “t” Xj)(y.i^j 0, A^ Aj. 


Divided by ai^ and sum those equations in (|6.11|) with indices (i, j) G X, we have 

n 

(6.12) fn n + N Afc = 0 . 

fc=i 

By incompressibility Afc = 0 and the global invariants asserted in ( |6.1CI|) follow. □ 


Two prototype examples are in order. In the 3D case we recover the global invariant 
(13), (Ai — A 2 )(A 2 — A 3 )(A 3 — Ai) = Const., corresponding to the sequence of pairs Xi = 
{(1, 2), (2, 3), (3,1)}. In the 4D case we have, in addition to incompressibility, the two 
global invariants. Hi := (Ai — A 2 )(A 3 — A4) = Consti and n2 := (Ai — A 3 )(A 2 — A4) = 
Const 2 , corresponding to the two X-sequences of indices, Xi = {( 1 , 2 ), (3,4)} and X 2 = 
{(1, 3), ( 2 ,4)}. Observe that a third global invariant corresponding to {(1,4), ( 1 , 3)}, is in 
fact generated by the difference of the hrst two, namely n3 = (Ai —A 4 )(A 2 — A3) = Hi — n2. 
Our next issue is therefore, a proper counting of these global invariants. 


6.2. On the number of global invariants. The proof of Lemma |6.2| makes clear the 
direct linkage between each X-sequence of indices, if,]) satisfying (| 6 . 8 |), and a global 
invariant formed by the corresponding product, n(jj)(Aj — Aj). Of course, not all the 
different X-sequences satisfying ( | 6 . 8 |) should be counted, since some of them lead to the 
same invariant products. We also need to remove any redundancy due to linear and 
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nonlinear dependence among these different invariant products. Thus, we inquire about 
the following 

Question. How many independent products, n(jj)gx(Aj — Xj) can be formed by X- 
sequences, i.e., sequences of indices {i,j) satisfying (|OD , 

n 

3N = N{n) e Z s.t. ^ ' 

(ij)ex k=i 

We know that the number of such independent invariant products together with the 
incompressibility constraint (|6.9|) does not exceed n, the number of independent global 
invariants of the restricted Euler (|6.1|) , and hence there are no more than n — 1 independent 
invariants of form ( |6.10|) . But the precise answer remains open, and in particular we are 
not clear whether all the global invariants of (| 6 . 1 | ) are necessarily the products formed in 
Lemma |6.2| . Below we provide a lower bound for our question, by the construction of 
such independent invariants. 

Let us begin by referring to the 4D example mentioned above. Starting with the first 
invariant. Hi = (Ai —A2)(A3 —A4), we derive a second independent invariant by exchanging, 
2 3, which leads to n2 = (Ai — A3)(A2 — A4). Other possible exchanges are redundant, 

say 2 ^ A yields the linearly dependent product n3 = (Ai — A4)(A2 — A3) = Hi — n2, and 
this is consistent with the fact that i = 3 , j = 4 playing a symmetric role in the original 
Hi-pair (A3 — A4). We conclude that while forming the linearly independent products, 
n(j j)gj(Aj — Aj), at most one ‘admissible’ exchange between different pairs of X-indices 
is permitted. Moreover, we should also exclude nonlinear dependence. For n = 8, for 
example, consider the four products. Hi = (Ai — A2)(A3 — A4)(A5 — A6)(A7 — As), n2 = 
(Ai — A3)(A2 — A4)(A5 — A6)(A7 — As), Hs = (Ai — A2)(A3 — A4)(A5 — A7)(A6 — As) and 
n4 = (Ai — A3)(A2 — A4)(A5 — A7)(A6 — As). The four invariant products are linearly 
independent — indeed, Q^fchlfc = 0 with A7 = Ag is reduced to a linear combination of 
the last two 4D independent pairs, asHg + a4n4 = 0 = 04 = 0 , and similarly, 

setting Ag = Ag yields ai = a2 = 0 . Nevertheless, they are redundant in view of their 
nonlinear dependence, n4 = n2 x Hg/Hi. 

Our construction of independent invariants in the general n-dimensional case proceeds 
as follows. We start, for even n = 2m, with the usual ordering Xi = (1, 2) (3,4)... (n — 
l,n). Making an admissible exchange between the first and second pairs yields the next 
independent invariant associated with X 2 = (1, 3) (2,4), (5, 6 )... (n — 1, n). Next, we make 
an admissible exchange between the second and third pairs, X 3 = (1, 2) (3, 5) (4, 6 )... (n — 
l,n), and so on. In this manner we proceed with one admissible exchange between each 
two consecutive pairs, leading to the m global invariants of the restricted Euler equations 

im 


Hi (Ai — A2)(A3 — A4) ■... ■ (A„_i — A„), 

Hfc := (Ai — A2) ■ ... ■ (A 2 fc -3 ~ A2fc-l)(A2fc-2 — X 2 k) ■ ■■■ ■ (An-l, A„), 
k = 2,3,... ,m. 

To verify that these m = | global invariants are independent, we note that by setting 
A 2 fc -3 = A 2 fc -2 we have H^- = | “ q’ 

Hfc ^ X’(ni,... , Hfcii,... , n„), k = 1,2,.. .m. 


(6.13) 


, which excludes the possible dependence 
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A similar procedure applies to the odd case, n = 2m +1. Starting with the usual ordering 
Xi = (1, 2), (2, 3),... (n, 1), we make an admissible exchange between the hrst and third 
pairs, I 2 = (1, 3)(2, 3)(2,4)..., the third and hfth pairs, (1, 2)(2, 3)(3, 5)(4, 5)(4, 6 )... 
and so on, leading to the m independent global invariants 


(6.14) 


Hi 

hlfc 


(Ai — A2)(A2 — A 3 ) ■... ■ (A„ — Ai), 

(Ai — A2)...(A2A:-3 ~ A2 fc-l)(A2 A:-2 ~ A2fc-l) (A2 A:-2 ” A2fc)...(A„ — Ai), 

fc = 2 ,3,... , m. 


We conclude with 


Lemma 6.3. [Global Invariants] The n-dimensional restricted Euler system ( \6.1\ ) subject 
to the incompressible initial data, ^ 11=1 -^*(0) = 0? admits the following + 1 global 
invariants in time: the incompressibility (\6.1C\) , — 0? additional n/2 

(- respectively, {n — l)/2) invariants specified in ( \6.13D for n even (and respectively in 
( \ 6 .i 4 ) for n odd). 


6.3. Behavior at the finite breakdown time. The rest of section is devoted to study 
the topology of the flow at the breakdown time based on the Lemma |6.2| . We start by 
noting that the level set of the integrals of the restricted flow 


(| 6 . 10 |) are not compact, 
and hence we have to perform singularity analysis to hgure out in which orthant the flow 
may diverge. The idea is to build local solutions around the singularities in order to 
study the blow up-rate and the location where the hnite-time blow-up actually occurs. 
The singularity analysis is a standard method to prove the integrability of ODEs. For 
readers’ convenience we sketch the main steps below, and refer to 0 and references 
therein for more details of this method. 

We assume a flow governed by the nonlinear ODE w' = f{w), diverges at a hnite time 
t*, and we then seek local solutions of the form 


W = UJT^ 




ajT 


Hi 




where t = t* — t, p & IR”, g G IN and Oj is a polynomial in log{t* — t) of degree Nj < j. 
There are three steps to determine the above series : (1) hnd the so called balance pair, 
{u},p), such that the dominant behavior, is an exact solution of some truncated 

system w' = f{w)] ( 2 ) computation of the resonances, which are given by the eigenvalues 

of the matrix —— diag{p); (3) the last step of the singularity analysis consists of 
Ending the explicit form for the different coefficients aj by inserting the full series in the 
original system, w' = f{w). 

Armed with the above algorithm, we proceed to carry out the singularity analysis for 
the restricted Euler equations. Let the dominant behavior of the A-system ( 0 ) assume 
the form 


Ai(t) 


~ (U.-T' 


■Pi 


Upon substitution into (|6.2|) one has 


i = 1, • ■ ■ ,n. 


-uJiPiT^^ ^ = 


n 

-E 

n ^ 

i=i 
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Equating the powers of r as r —>• 0 we find, pi = —1, and the ends satisfy the equation 




1 

n 



There is a n-parameter family of such cd’s 


n- 2 ’ ’n- 2 ’ ’n -2 


k = 1, ■ ■ ■ ,n. 


Due to the symmetry of the equation, the flow may diverge in n out of 2”' orthants. More 
precisely, we have 


Lemma 6.4. The onlyn stable solutions of the spectral dynamics ^6.^) , A = (Ai, • • • , A„), 
associated with restricted Euler equations (H are explicitly given by 

A = (n-2)a(a;) ^ a('=)(0 ) = = 

n — 2 — a[x)t 

with arbitrary a{x) < 0 . 


To sum up, we state the following 


Theorem 6.5. Consider the restricted Euler dynamics l\ 6 . 3 i ) with initial data (Ai(0), ■ ■ ■ , A„ 
The level set of global invariants given by 

nijGx(Ai - Xj) = Const 


are not compact. The general solution may break down at finite time in one of the n 
orthants {+, +, • • ■ ,—,•••+,+} along the kth separatrix 


(I,-- 


whenever a(x) > 0 . 


,1 


n, 


, 1 ) 


a(x) 


n — 2 — a(x)t’ 


k = 1, - ■ ■ ,n 


Remark 6 . 6 . Other possible variants of the local restricted Euler equations can be written 
in the form 

J^M + e{M^ - trM’^jnln^n) = 0 

with d G (0, cxo). This equation becomes anisotropic, but the local topology of the solution 
remains the same as in the isotropic model ( |6.1D below. Indeed a hyperbolic scaling, 
{t,x) —>■ {9t,9x), leads to the isotropic model corresponding to 6 ^ = 1 . 


7. Restricted Euler-Poisson dynamics 


We begin by introducing the so called restricted Euler-Poisson equations. As argued 
in we retain the local part of the nonlocal term kR[p] in the Lagrangian form of 
the Euler-Poisson (EP) system (|2.17|) , (|2.18|) to obtain a restricted Euler-Poisson system 
( CT ), (iH, i.e.. 


dfM -|- M ■ VAf -|- — —pinxni 

n 

dtp -I- u ■ Vp -I- ptrM = 0. 
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If we let Aj(a;, t) denote the eigenvalues of velocity gradient tensor Vn, then by the spectral 


dynamics Lemma p.l|, the eigenvalues and the density p are coupled through 


(7.1) 


(7.2) 


dtp + u - Vp + p ^ Xj = 0 , 
i=i 

dtXi + M ■ VAj + ^ 


i = 1, 


,n. 


n 


This is a closed system governing the restricted Euler-Poisson equations, which serves as 
a simple approximation for the evolution of the full Euler-Poisson system (|2.14|) - (p.l6|) . 

In this section we use the spectral dynamics of the restricted Euler-Poisson equations 
to show two main points: 

1) . The global existence of the smooth solutions for a large class of 2D initial conhgurations- 
consult Theorem o below. 

2) . The hnite time blowup of the n-dimensional solutions subject to another class of 
initial data outlined in Theorem |7.4| below. As a consequence of 1) and 2), it follows 


that the 2D restricted Euler-Poisson equations admit a critical threshold which distin¬ 
guishes between initial conhgurations leading to hnite time breakdown and global smooth 
solutions. A detailed study of this 2-dimensional critical threshold phenomena in this 
context is provided in [^. This complements the study of critical threshold phenomena 
for isotropic conhgurations in the general (global) Euler-Poisson equations presented in 

0 . 

We start with the global regularity of 2-D restricted Euler-Poisson solutions. By well 
known arguments, the global regularity follows from local existence complemented by 
a boot-strap argument based on the apriori estimate of ||Vn||Loo. For the 2D restricted 
Euler-Poisson model, the velocity gradient tensor Vu is completely controlled by its eigen¬ 
values, Ai,f = 1,2, consult for a detailed statement of this argument. With this in 
mind, it is left to obtain apriori uniform bound of Aj’s yielding a sufficient condition for 
the global existence of smooth solutions for the restricted Euler-Poisson model. 

Theorem 7.1 (Global existence). The solutions of the 2-D restricted Euler-Poisson equa¬ 
tions &([Z3) remains smooth for all time t > 0 if both Aj(0), i = 1,2 are complex, i.e., 
Im{Xi{a,0)) ^ 0. 

Proof. In the 2-D case the density equation o becomes 

p' + P(-^i + -^ 2 ) = 0, ' := dt u ■ V^. 

/^From (7^) it follows that the evolution of the divergence d = Ai -|- A 2 , is governed by 

d' -\- d'^ — 2 X 1 X 2 = kp, 
and the evolution of A = A 1 A 2 is given by 


(7.3) 

Introduce the ‘indicator’ function. 


k 

A' -\- dA = -^pd- 


r(t) := exp( / d{x{a,T),T)dT), 
Jo 
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then the density equation gives 

p{x, t) = po{a)/T{a, t), t > 0. 

Noting that 

r' = dr, r" = + d')r, 

we then have 

(7.4) r" - 2Ar = kpo. 

Substitution of d = T'/F and p = po/T into the A equation (|773|) it follows that 


(TAy = ^[l„T]’. 


Integration once gives 


kpo. 


TA — —^InT + Aq, Aq (AiA2)(a, 0), 


which when inserted into ( [f.4|) yields 

r" = kpolnV + 2Ao + kpo - 

The integral energy becomes 

[r^]^ = dg + 2(2Ao + kpo){T — 1) + 2fcpo J In^d^ 

= dg + 4Ag(r — 1) + 2 kpoTlnT . 

Assume that the solution break down at a hnite time t*, i.e., r{t*) = 0, then at this time 
one has 

[r']^ = dg — 4Ag = (Aio — A2 o)^- 

Therefore hnite time breakdown can not occur if Ai(a, 0) is complex. □ 


Remark 7.2. The above sufficient condition is satished, for example, by the initial ve¬ 
locity with large enough vorticity oj := — Vy, associated with the scaled velocity 

{uo{l3x,y),Vo{x,py)) with (Mg(j, ■),T)) such that |a;o| ~ » |do|, implying dm (A ig) 7 ^ 

0 . 


Remark 7.3. What happens with a possible blow-up if both Ajo are real? Let t* be a hnite 
blow-up time. If follows that the blow-up rate is necessarily of the form {t* —i.e., 

d{t) ~ ^— as ft t*. 

f t 

This follows from a simple analysis on the following relations 

0 > r'(f*) = Aio — A 20 = lim d(f)r(f), r(f) = exp( [ d{T)dT). 

do 
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The loss of smoothness of the velocity held is closely related with the intricate problem 
of weak convergence in the absence of the strong convergence. The open question in this 
context is how the nonlocal term affects the topology of the how. 

To gain further insight on the question of global regularity vs. hnite time breakdown, 
we continue with the u-dimensional restricted Euler-Poisson dynamics 
before, we subtract two consecutive eigenvalue equations in (|7.2|) to obtain 

[ln{Xi — Aj)]' + Aj + Xj = 0, for i ^ j. 

Summation over (i, j) G X, with X dehned in (| 6 . 8 |) , gives 

n 

[Z?7,(nj(j)£j(Aj — Aj)]^ + N y Afc = 0, 

k=l 

Combined with the density equation, [Inp]' + Yyk=i = 0, this yields the following global 
invariants 

( 7 . 5 ) Ep,,).xjA, - A,) ^ 

The level set of the above invariants is not compact and the hnite time singularity can 
not be ruled out, and in fact, noncompactness implies that certain portion of the phase 
space must lead to hnite time breakdown. 

In order to perform a singularity analysis similar to the one provided in §P, we consider 
a truncated system 

n 

p = -p^K 

i=l 

K = -^i * = 1 , • • • , n. 


(7.1)-(7.2). As 


Finding its dominant solution of the form 


(Pj ‘ ) Xn) ^ (^0) ^1) ■ ■ ■ 5 

with p = {po,pi, ■ ■ ■ ,Pn) and t = t* — t leads to 


-UJqPqTP^ ^ 


-OOiPiT^^ ^ 


-a;or^° 


—ojJt 


“ipi 


This gives the balance (a;,p) with 

a; = (1, -1, ■ • • , -1), p = (-U, -1, • • • , -1), 

where cuo > 0 is chosen so that it is consistent with the positivity of the density. 

Therefore, there exists a general Puiseux-series solution |]^ based on the above balance 
pair, and the blow-up may occur on the orthant —, • • ■ , —This, combined with the 
noncompact integrals derived in (|7.5|) , shows that the solution must exhibit hnite-time 
blow-up in the above orthant. 

To summarize, we state the following. 
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Theorem 7.4 (Global invariants for n > 2). Consider the restricted Euler-Poisson dy¬ 
namics with real initial data (po, Ai(0), • ■ ■ , A„(0)). Then, there exist global 

invariants 


(7.6) 


bl(i,i)Gx(Ai Xj) 


P 


N 


Const., N 


1, n even, 

2, n odd. 


The general solution may break down at finite time in the orthant {+, —, • • • , —}■ 


Two particular cases are worth mentioning. In the 2-D case we have one global invariant 
(A 2 — Ai)/p = Const., while the global invariant in the 3-D case, corresponding to X = 
{(1,2), (2, 3), (3,1)}, is given by 


(Ai — A2)(A2 — A3)(A3 — Ai) 


Const. 


8. Concluding Remarks 


This work provides a general framework for several variants of the restricted Euler 

m 


dynamics in multi-dimensional case, extending the previous study initiated in 


The 


main tool in this paper is the spectral dynamics analysis. We should point out that this 
analysis enables us to derive global invariants which are otherwise difficult to detect— 
one such example was used with the viscous dusty medium model in §5. In particular, 
we obtain a family of global spectral invariants, interesting for their own sake, for both 
restricted Euler equations (|6.9|) - (|6.1(]|) and the restricted Euler-Poisson equations (|7.6|) . 

Noncompactness of the level set of these global invariants implies the finite time break¬ 
down for a class of initial configurations, for which the local topology of the restricted flow 
is analyzed. This was demonstrated in Theorem p.5| in the context of the restricted Euler 
equations. The finite time breakdown in this restricted model does not necessarily bear 
on the full, non-restricted Euler equations. On the other extreme we have the possible 
scenario of a global existence of smooth solutions for restricted models such as restricted 
Euler-Poisson equations, for which we have the global existence once a critical threshold 
condition is met. Here, we believe, the global existence does carry over to the question 
of global existence for the full non-restricted Euler-Poisson equations. In particular, in §7 
we have shown the existence of a critical threshold for the 2D restricted model, which in 
turn suggests the critical threshold phenomena for the full 2D Euler-Poisson equations. 

We close this section with following comments. 


Remark 8.1. Suggesting other nonlocal restricted models. To gain further insight on the 
fine structure of the flow we propose the following restricted nonlocal models for both 
Euler-Poisson equations and the incompressible Euler-equations, the analysis of which 
will appear elsewhere. 

• Euler-Poisson equations 

We take the diagonal part of the right side of the M-equation in the Euler-Poisson 
dynamics (|2.17| ), (|2.18|) to obtain 


( 8 . 1 ) 

( 8 . 2 ) 


dtp -I- u ■ Vp -I- ptrM = 0, 

dtM + u ■ VM + = k{RiRj{p)6ij). 
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• Restricted Euler-dynamics 

A restricted nonlocal Euler dynamics 

(8.3) dtM + u ■ VM + 

We note that of course dttrM = 0 and the incompressibility is still invariant. 


Remark 8.2. Is the spectral dynamics sufficient? We are aware that the spectral dynamics 
does not tell the whole story for general fluid flows. The following example of a Burgers 
shear-layer pm demonstrates this point. Here the simplest solutions of the inviscid 
Euler equations are the Burgers shear-layer solutions with the velocity held given by 

u = ih{x 2 , t), 0, 0) + (0, -7x2, 7 x 3 ). 


The velocity gradient tensor, 

/ 0 0 \ 

Vu = 10 — 70 , 

V 0 0 7 / 

has eigenvalues {— 7 , 0 , 7 } which rehect strain effects, but otherwise are independent of the 
arbitrary shear-layer effect h{x 2 ,t)- Thus, the eigenvalues can not capture the complete 
behavior of this h—dependent how. 


Remark 8.3. A main issue in this context is how the restricted Euler-type dynamics re¬ 
lated to real hows and at what scale of motion it might apply. In the Navier-Stokes 
equations, for instance, the nonlocal term should not be ignored at both large and small 
scales. At large scales the pressure-driven eddy intersections are important and at small 
scales the velocity gradients are limited by viscous dihusion. We refer the reader to [|I 2 


for a detailed discussion on this issue. Another interesting issue left for future research is 
the recovery of the gradient velocity tensor from the known spectral dynamics. 


9. Appendix. Trace dynamics for the restricted Euler-equations 

This appendix is devoted to an alternative formulation of the spectral dynamics in 
terms of the traces of M^, k = 1, ■ ■ ■ ,n, where M solves the restricted Euler equation 


(9.1) 


+ tr = i(rAfV„ 

at n 


This is motivated by the trace dynamics originally studied in [Q for n = 3. 

Here we seek an extension for the general n-dimensional setting, which is summarized 
in the following 

Lemma 9.1. Consider the n-dimensional restricted Euler system B) subject to the 
incompressibility condition mi := trM = 0. Then the traces m^ := trM^ for fc = 2, • • • , n 


satisfy a closed dynamical system, see B)-B) with ( \9.(\ ) below, which governs the local 
topology of the restricted flow. 
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Proof. Based on the equation ( |9.1|) the transport equations for higher products of M can 
be written as 

+ 2M^ = -MtrM\ 
at n 

+ 3M^ = -MhrM^, 
at n 


at 


Taking the trace of the above equations and using m 2 = with mi 


(9.2) 

(9.3) 


d 

—1712 + 2m3 = 0, 

dt 

d 3 

dt n 


2 

2 ) 


0 leads to 


(9.4) —m„ + nm„+i = m„_im2. 

To close the system, it remains to express rrin+i in terms of (mi, • ■ ■ ,mn). To this end 
we utilize the Cayley-Hamilton theorem 

(9.5) + giM"-' + ■ ■ ■ g„_iM + qj = 0, 
expressed in terms of the characteristic coefficients 

gi = -mi = 0 , q2 = ~m2, ga =-"13/3, ^4 =-^^ 4/4 + m^/S, ■■■. 

Note that the g’s can be expressed in terms of (mi, ■ ■ ■ , m„). (The procedure for comput¬ 
ing these coefficients is given at the end of this appendix). Using the Cayley-Hamilton 
relation ( |9.5|) one may reduce m„+i in (|9.4|) to lower-order products. In fact, tr{M x ( |9.5| )) 
gives 

(9.6) m„+i -F q2mn-i H-h qn-im2 = 0. 

Substitution into (|9.4| ) yields the closed system we sought for. □ 

We now turn to consider two examples which demonstrate the above procedure. 


Example 1. (3-dimensional case n = 3, see 
In the three dimensional case one has 


hence 


qi =0, q 2 = - 2 ^ 2 , qs = det{M) 


o 1 1 

M'^ - m2M -m 3 = 0 . 

2 3 





This gives 
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Thus a closed system is obtained, 

(9.7) 

(9.8) 


d 

—1712 + ‘^rns 

dt 

d In 

—rri'i H —rrin 

dt ^ 2 ^ 


0 , 

0 . 


The invariant of 6 m| = Const., could be easily obtained. We consider the phase plane 
(m 2 , m 3 ), except for the separatrix 6 m| = m|, all other solutions would not approach the 
origin and have the hnite time breakdown, see Figure E3 



Figure 9.1. 3-D Trace-dynamics in Restricted Euler Equations 


Example 2. (4-dimensional case) 

In the four dimensional case one has 


1 1 7714 mi 

qi = 0 , q2 = -^m2, qs = Qa = —^ + —• 


Hence 




-m2M^ -msM 


m 4 m^ 


2 3 4 

Multiplying by M and taking the trace we have 




1 


1 


ms = -m2m3 -I— m^m2 = -m2m^. 
2 3 6 


Therefore the resulting closed system becomes 


(9.9) 

(9.10) 

(9.11) 


d 

—m 2 + 2 m 3 = 0 , 

dt 

d ^2 

—m 3 -F 3 m 4 = -m 2 . 


d 


7 


—m 4 = —m3m2. 

dt 3 







28 


H. Liu and E. Tadmor 


This system is still integrable with the following two invariants 

3C 

3ml = ml H— ^"^2 + C2, 12m4 = 7ml + Ci, 

where Ci, C 2 are constants integrals of the flow. 


Remark 9.2. Note that when Ci = 0, the projection of the trajectory on the m 2 —m 3 
plane has the same topology as that in the 3-D case. See Figure p.2| for the vector field 
in (m 2 ,m 3 ,m 4 ) space. 


4-D Trace dynamics 



Figure 9.2. 4-D Trace-dynamics in Restricted Euler Equations 


To gain further insight on the formation of the singularity in this case, we try the 
dominant solution of the form with t = t* —t, for the truncated system 

d 

—m 2 = -2m3, 
dt 

d ^2 

= 4"^- 


A simple computation gives 


d 


7 


—mi = --m3m2. 

dt 3 


p =(-2,-3,-4), a = (4,-4,y), 
which shows that the flow may diverge in the orthant 


Remark 9.3. The above examples demonstrate the difficulty in deriving the global invari¬ 
ants for arbitrary n > 3 equations, without the insight provided by the spectral dynamics. 


We now conclude this appendix by presenting a procedure of computing the coefficients 
in the characteristic polynomial for a given matrix. 
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Lemma 9.4. Let A be a square matrix of order n, its characteristic polynomial reads 

n 

det{XI - A) = 

k=0 

Then qj = tr{AyA)), where A'^(A) is the tensor product of A. 

Proof. Noting that for e = — 

n 

det{I + eA) = y2i-^)"(lke\ 

k=0 

On the other hand, if |e| small, then 

trlog{I + eA) = 

converges. These two relations when combined with the identity 

det{I + eA) = exp {trlog{I + eA)) 


yield 


= exp 




fc =0 


vfc=l 


k 


Eqnating the same powers of e on both sides gives 

% = 1, 

qi = -ai, 

o 

02 


Uj2 U/i 

03 CI1CI2 

® = “T + ~ 

O4 CI1CI3 

94 = - j + — 


3!’ 


+ — 
8 


0^02 


Ui 

+ 4! ■■■’ 


where Ok = trA^. This procednre gives the expression of each qj in terms of for 
= 1, ■ • ■ , n. □ 
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